SLAC-PUB-10680 
hep-ph/0408344 

August 30, 2004 



Loop Corrections to Heavy-to-Light Form Factors and 
; Evanescent Operators in SCET 

o 
o 

< 



Thomas Becher and Richard J. Hill 



Stanford Linear Accelerator Center, Stanford University 
Stanford, CA 94309, U.S.A. 

m 
oo . 

2 ' Abstract 
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One-loop matching corrections are calculated for Soft-Collinear Effective Theory (SCET) 
operators relevant to the analysis of heavy-to-light meson form factors at large recoil. 
The numerical impact of radiative corrections on form factor predictions is assessed. 
Evanescent operators in the effective theory are studied and it is shown that even in 
' problems of the Sudakov type, these operators can be renormalized to have vanishing 

matrix elements. 
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1 Introduction 



The methods that were used to establish factorization theorems for high-energy QCD processes 
can also be applied to study decays of 5-mesons into light hadrons. In this case the hard scale 
is set by the mass of the heavy 6-quark and factorization theorems for inclusive [1, 2, 3] as 
well as exclusive [4] 5-decays to light hadrons arise in the heavy-quark limit. Soft-Collinear 
Effective Theory (SCET) is the effective field theory arising in this limit. It describes the 
heavy-quark expansion for these decays, and permits the study of their factorization and 
renormalization properties [5, 6, 7, 8, 9]. 

Semileptonic decays, such as B — > niu, are the simplest exclusive heavy-to-light meson 
processes. These decays are described by weak-interaction form factors, which at large recoil 
energies of the light meson, E~m B /2> Aq CD , take the form [10], 

F t B ^ M (E) = C l (E)( M (E)+ / — 0bH / dufuMuWE^u), (1) 

Jo u Jo 

up to corrections suppressed by Agco/m^. Here M denotes the final state (pseudoscalar 
or vector) meson. The process-dependent Wilson coefficient functions Ci(E) and Ti(E,u,u) 
are calculable in perturbation theory and are the subject of the present paper. The wave 
functions, or more precisely, light cone distribution amplitudes (LCDAs) 4>b{w) and <Pm{u)-, 
and the functions (m(E), are process- independent nonperturbative hadronic parameters. The 
factorization theorem (1) has recently been established in the context of SCET [11, 12, 13]. 

At large recoil, the form factors Ef^ M are in fact subleading quantities, requiring power- 
suppressed interactions to mediate the B — ^ M transition. Since the light degrees of freedom 
inside the B- meson carry soft momenta, p M ~ Aqcd, whereas the degrees of freedom inside 
the final-state meson typically share the large energy of the meson, the transition can only 
occur if: (i), the partons are in an atypical "endpoint" configuration, allowing the soft B- 
meson constituents to be absorbed into the energetic light meson; or (ii), a large momentum 
is exchanged between the active quark and the light spectator degrees of freedom. Possibility 
(i), the so-called soft-overlap mechanism, is described by the non-factorizable part (m(E) of 
the form factors, and in the effective theory is given by the matrix element of the leading- 
order heavy-to-light quark current. This current involves only the two-component spinor 
fields of SCET and Heavy-Quark Effective Theory (HQET) [14], and mediates the transition 
of the heavy quark inside the 5-meson into an energetic light quark component of the light 
meson. The simple spin structure of the two-component fields results in "large-energy" spin 
symmetries, such that a single function Cm(E) describes the soft-overlap contribution to all 
form factors involving the same final-state meson M [15]. Possibility (ii), the hard-scattering 
mechanism, results from subleading effective theory current operators involving an additional 
gluon field transferring large momentum to the spectator. Their matrix elements factorize 
and correspond to the second term in equation (1). 

The form factors receive contributions associated with different momentum scales: the hard 
scale /i 2 ~ m 2 , an intermediate scale fx 2 ~ m fe AQ CD and the soft scale /i 2 ~ Aq CD . In order to 
perform a renormalizat ion-group (RG) resummation of large perturbative logarithms involving 
ratios of these scales, a two-step matching is required, first from QCD onto an intermediate 
theory, denoted SCETi, and then from SCETi onto the final low-energy theory, denoted 
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Figure 1: Two step matching from QCD to SCETi to SCETn. The double line denotes 
the heavy quark; dashed lines denote hard-collinear fields in SCETi and collinear fields in 
SCET n . The symbol "<g>" indicates the convolution of the Wilson coefficient functions with 
the non-local operators. 



SCETn; see Figure 1. For the leading-power currents defining the soft-overlap contributions, 
the matching from QCD onto SCETi has been performed through one-loop order in [5], and 
the subsequent matching onto SCETn is given in [13]. Resummation of large logarithms in 
the coefficient functions Ci(E) results in a universal RG factor [5, 13], the same for all form 
factors. However, the numerical value of this universal factor is presently of limited utility, 
since the functions (m{E) are not known at any fixed renormalization scale. On the other 
hand, some predictions exist for the hadronic wavefunctions appearing in the second term 
of (1), and we may study the numerical impact of renormalization-group evolution down to 
hadronic scales where these predictions are reliable. In a recent paper in collaboration with Lee 
and Neubert [16], we resummed all leading single and double logarithms in the hard scattering 
kernels Ti(E, u, u). At leading order in RG- improved perturbation theory, the resummed result 
is obtained using tree-level matching. However, since a s is sizeable at the relevant scales, we 
expect one-loop matching corrections to also be important. In [16] the size of these matching 
corrections was investigated using the scalar current as an example, and they were found to 
be comparable to the resummation effects. In the present paper, we compute the complete 
one-loop matching results for the hard-scattering term, for both matching steps. 

The outline of the paper is as follows. In Section 2 we give the matching coefficients 
for QCD vector and tensor currents onto the subleading SCETi operators relevant to the 
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hard-scattering form factor contributions. The coefficients for the scalar current have been 
calculated previously [16, 17], and can in fact be derived from the vector current results [16]. 
The vector and tensor coefficients have also been calculated by Beneke, Kiyo and Yang [17] and 
we find full agreement with their results; in Section 2 we list only those coefficients necessary 
for B — > M form factors, in the "primed" operator basis of [16]. The Wilson coefficients, 
also called jet functions, that arise in the second matching step, from SCETj onto SCET n , 
are considered in Section 3. We have already stated the results relevant to the form factors 
in [16]. Here we present the general expressions from which those results were derived, and 
provide details of the calculation. 

For dimension d = 4 — 2e in dimensional regularization, so-called evanescent operators [18, 
19, 20] arise in the matching of QCD onto SCET. The matrix elements of these operators 
vanish at tree-level for d — 4, but in a generic renormalization scheme they are nonzero at 
higher order. To avoid introducing new hadronic functions at each order in perturbation theory 
to parameterize the matrix elements of the evanescent operators, it is important to show that 
a renormalization scheme exists for which these operators vanish exactly in four dimensions. 
The existence of such a scheme is related to the technical observation that in loop Feynman 
diagrams involving an evanescent operator, a contribution proportional to a physical operator 
contains a factor e from Dirac algebra, which cancels a 1/e divergence arising from integration 
over loop momenta; the resulting finite contribution is local, and so may be subtracted by 
a suitable counterterm. It can also be shown that in this scheme no mixing occurs of the 
physical into the evanescent operators, and thus the matching coefficients onto renormalized 
evanescent operators are irrelevant to the calculation of physical matrix elements. In processes 
involving Sudakov double logarithms, the existence of 1/e 2 divergences in loop integrations 
poses a potential obstacle to these standard arguments. In Section 4 we show that despite 
this complication, renormalization schemes exist in which the matrix elements of renormalized 
evanescent operators vanish. This section also discusses a related issue, the choice of evanescent 
operators. Different choices correspond to different renormalization schemes for the physical 
operators. For the four-quark SCETn operators describing the hard-scattering form factor 
contributions, we isolate the particular operator basis which corresponds to the MS scheme 
after Fierz transformation. Upon taking heavy-light meson matrix elements, this ensures that 
the LCDAs appearing in (1) are renormalized in the MS scheme. 

The phenomenological importance of one-loop matching corrections to heavy-to-light form 
factors is investigated in Section 5. Hard-scale corrections, corresponding to QCD — > SCETj 
matching, are found to be of order 10 — 20%. These corrections determine the size of violations 
to the large-energy spin-symmetry relations obeyed by the soft-overlap term, and also the size 
of non-universal corrections to the hard-scattering term. We notice that in the soft-overlap 
case, the radiative corrections for different form factors are remarkably similar, and therefore 
have little effect on the symmetry relations. The numerical impact of one-loop contributions 
to the jet functions for SCETi — > SCETn matching is more difficult to estimate, owing to the 
lack of precise knowledge on the _B-meson wavefunction. However, when hard-scale QCD — > 
SCETj matching corrections are ignored, the hard-scattering contributions are described by 
the universal functions Hm{E) introduced in [16]; radiative corrections to the jet functions are 
therefore not relevant to the conclusions which can be drawn based solely on this universality. 
The radiative corrections are necessary in order to relate the functions H M {E) to models for 
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the 5-meson wavefunction, and we investigate the size of the corrections for two such models, 
finding that they increase H M (E) by ~ 20 — 30%. Section 6 provides a discussion and our 
conclusions. 

2 One-loop matching from QCD onto SCETj 

In this section, we deal with SCET^ which describes the interaction of energetic light particles 
of virtuality p 2 ~ m b Aq CD with a static heavy quark. The soft sector of SCET^ is described 
by HQET, with a heavy quark of velocity v. Different components of hard-collinear momenta 
and fields scale with different powers of the expansion parameter A 1 / 2 = (Aqcd /^b) 1 ^ 2 - The 
large momentum components are isolated to obtain a definite power counting by working with 
light-like reference vectors, n M in the direction of the jet of outgoing collinear particles, and 
n M a complementary vector satisfying n ■ n = 2: 

P fl = n-p— +n-p— +p^_. (2) 

The components (n-p, n-p,pj_) of a hard-collinear momentum are defined to scale as (A, 1, A 1 / 2 ), 
while for a soft momentum the components scale as (A, A, A). We choose n M such that v± = 0, 
implying n M = (2v^ — n^/n-v)/n-v. The canonical choices are t> M = (1, 0, 0, 0), tV 1 = (1, 0, 0, 1) 
and n M = (1, 0, 0, —1), corresponding to the energetic jet in the z-direction as viewed from the 
rest-frame of the heavy meson. Throughout this paper we use the notation and conventions 
of [16], to which we refer the reader for details and references. 

We now consider the effective theory representation of the scalar, vector and tensor flavor- 
changing currents, 

S = qb, V^ = q^b, = (-i)qa^b = q^Y ] b- (3) 

We use square brackets around indices to denote antisymmetrization. The QCD scalar and 
tensor currents require renormalization. We work in dimensional regularization for d = 4 — 2e 
dimensions, and define the QCD as well as the effective theory operators in the modified min- 
imal subtraction (MS) scheme. 1 We use the naive dimensional regularization (NDR) scheme 
for 75 so that the pseudoscalar and pseudovector currents renormalize identically to the scalar 
and vector currents, and our results for the matching are independent of the chirality of the 
light quark q in the current operator. The construction of the effective theory currents is 
discussed in [16]. For the case of the scalar current, which we take at position x = 0, only two 
effective theory operators are needed through the first subleading order: 

J|( S , x = 0) = X hc (sn) (l - ^= f) h(0) , 

V in - d 1 J 

J#(s, r,x = 0) = X hc {sn) 4 hc± {rn) h(0) . (4) 

1 In Section 4, we consider the modifications to this scheme necessary to ensure the vanishing of renormalizcd 
evanescent operators. 
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The field h is the heavy-quark field, X^ c = W^hc is the hard-collinear quark field multiplied 
by a hard-collinear Wilson line Wh c in the n-direction, and A^ c± = W^ c [{id^_ + ^ c j_)W/i C ], 
where Ah c is the hard-collinear gluon field. Note that the above operators are non-local: 
the hard-collinear fields live at different points on the light-ray in the n-direction. This is 
a consequence of the fact that n ■ d derivatives on hard-collinear fields count as quantities 
of order one, so that operators involving arbitrary functions of such derivatives may appear 
at a given order in the power counting. The Wilson coefficients of these operators are then 
functions of the position arguments and the representation of the scalar current at x = in 
the effective theory reads 

q b ^ J ds ci( s )ji(s) + ^ J drdsC B s {s,r)Jl(s,r) + ... 

= Cj(E) J#(0) + ± J duC§(E,u) Ji(u) + . . . , (5) 

where the factor 1/2E has been inserted into the second term for convenience. We use a tilde 
to denote the coefficients in position space. The momentum space coefficient functions in the 
second line are defined as 

Ci(E) = J dse^ p Cj(s), 

C§(E,u) = J drdse l(us+ar)n - p Ci(s,r). (6) 

Here P = P out — P m is the total hard-collinear momentum of external states, and E = 
(n • v)(n • P)/2. The variable u G [0, 1] denotes the fraction of the energy E carried by the 
quark field Xh c and u — 1 — u the fraction carried by A^ cL . 

The fact that the subleading term in Jg involving a perpendicular partial derivative has 
the same Wilson coefficient as the leading order term is a consequence of reparametrization 
invariance, the fact that QCD is independent of the reference vectors and used in the 
construction of the effective theory [7, 21, 16]. Due to the presence of the reference vectors 
and v^, the vector and tensor currents are represented by several effective theory operators, 
even at leading order. However the general structure of the subleading currents is the same 
as in the scalar case: A-type operators involve perpendicular partial derivatives and have 
the same Wilson coefficient as the leading order operators; the remaining, B-type operators 
involve the insertion of an additional perpendicular gluon field A^ c± . The complete basis of 
SCETj current operators can be found e.g. in [16]. In that reference two different sets of basis 
operators are used, Jf and .if . Here we use the £?-type operators in the primed basis, which 
are more convenient in that they renormalize multiplicatively and match in a simple way onto 
SCETn- Suppressing the "/ic" subscript on hard-collinear fields, the vector current operators 
in this basis are 

= X( S n) 4 ± (rn) 7^ h(0) , Jyf = X(sn) 4 ± (rn) v» h(0) , 

= %{sn) 4±{rn) — h(0) , = X(sn) 7^ 4±(rn) h(0) , (7) 
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Figure 2: One-loop QCD diagrams contributing to the matching calculation for the subleading 
scalar current. The double line denotes the heavy quark. 
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(8) 



One-loop matching corrections to the leading order, A-type effective theory currents were 
evaluated in [5] and the results were confirmed by the findings of [17]. Here, we focus on 
the Wilson coefficients of the 5-type operators which can be obtained by evaluating the 
QCD diagrams shown in Figure 2. By performing the matching on-shell, one can avoid the 
calculation of effective theory loop diagrams, as discussed in [16], where we presented the 
one- loop matching for the scalar current. Only six of the above currents match onto SCETn 
four-quark operators relevant for the heavy-to-light form factors and we now give the results 
for those six Wilson coefficients. Introducing x = 2E/m b , we find for the case of the vector 
current, 
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The three coefficients needed for the tensor case are 
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We have checked that our results for all four vector currents, and seven tensor currents (in- 
cluding those not listed here) are in complete agreement with previous results of Beneke, Kiyo 
and Yang [17] once we translate to their operator basis, using the relations 
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for the vector coefficients, and for the tensor coefficients: 
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3 One-loop matching from SCETj onto SCETjj 



We now map the -B-type SCETi currents onto SCET n four-quark operators, as represented 
in the right-hand side of Figure 1. We recall that SCETn is the low-energy effective theory 
encompassing collinear modes with momentum scaling in ■ p,n ■ p,p±) ~ (A 2 , 1, A), and soft 
modes with momentum scaling (A, A, A). In general, soft-collinear messenger modes, with 
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Figure 3: SCET^ graphs contributing to the matching of the currents Jf (crossed circle) onto 
the SCETn four-quark operators Oj. Full lines denote soft fields, dashed lines hard-collinear 
fields. Diagrams with soft gluons or scaleless loops are not shown. 



momentum scaling (A 2 , A, A 3 / 2 ) are also required, even at leading order, and may communicate 
between the soft and collinear sectors [22]. However, in the particular case of the color-singlet 
four-quark operators under consideration, the sum over all such contributions can be shown 
to vanish; this is equivalent to the statement that the soft and collinear sectors factorize for 
the corresponding matrix elements [23, 13]. The resulting factorizable contributions to the 
form factors take the form of the second term in (1). The function Ti(E,u,u) appearing 
in this formula is a product (in general, a convolution) of the SCETi coefficients Cf from 
the first matching step, considered in Section 2, and the Wilson coefficients occurring in the 
second matching step, called jet-functions. The fact that the heavy quark and the collinear 
quarks are described by two component spinors restricts the allowed form of possible SCETn 
operators, and leads to relations among the jet-functions [16]. For the operators relevant 
to 5-decay form factors, only two jet functions are relevant; these are denoted by J\\ and 
J7i, and describe decays into pseudoscalar or longitudinally-polarized vector mesons, and into 
perpendicularly-polarized vector mesons, respectively. For completeness, and to illustrate the 
effects of evanescent operators, we will present in this section complete one-loop results for all 
jet functions. 

The .B-type currents in SCETi may be written, 

J r (s, r) = X hc (sn) A hc±ll (rn) h(0) , (13) 

where the Dirac structure depends on the specific current under consideration; e.g. = 7^ 
for the scalar current. For the coefficient function of the operator J r (s,r), we write 

C r (u) = fdrds e ^ us+ar > p C r (s, r) . (14) 
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The four-quark SCET n operators depicted in Figure 1 take the form, 

O k r l ;r(s,t) = X k c (sn)T^K l (0) QT(tn)T^X n c (0) . (15) 

The superscripts k, I, m, n are color indices, which will be suppressed in the following. Since we 
are interested in matrix elements for color-singlet initial and final state hadrons, our primary 
concern is the color-singlet contribution, proportional to the contraction of (15) with S kn S ml . 
Here X c = W^ c is the restriction of X hc to collinear momentum modes. 'K = S^h and 
Q s = S^q s are soft fields multiplied by a soft Wilson line S in the n-direction. The precise 
definitions of these fields and the SCETn effective Lagrangian is given in [22]. 

Unlike in the first matching step, from QCD onto SCETi, the matching of SCETi onto 
SCETn generates an infinite tower of SCETn operators. From the Feynman rules of SCETi 
it follows that must take the form 

Fg£ — 7-L/i; 7_L/ii/i2M3 5 7X^1/12/^3/^4^5' • • • • (16) 

We denote totally antisymmetrized combinations of Dirac matrices by y^ 1 "-^™ = 7^ . . .7^. 
A complete basis is obtained by considering the most general T^ h with the appropriate number 
of Lorentz indices, constructed from and 7^ (recall that we choose the heavy-quark velocity 
such that = 0). At one-loop order, three structures are possible: 

® 7±m , i = 1 , 

i>m ® = { ^ ±p r p ± ® 7±M , 1 = 2, (17) 
7^r^ <g> ^ ±flup , % = 3 . 

Operators for which involves antisymmetric products of five and more Dirac matrices 
appear only at two and more loops. In four dimensions, antisymmetric products of more 
than two perpendicular Dirac matrices vanish. Operators such as Or,3(s,t) containing these 
Dirac structures appear at one-loop and beyond in dimensional regularization, and are called 
evanescent [18, 19, 20]. The choice of evanescent operators is not unique, and affects the 
coefficients of the physical operators. We will discuss this issue, and the related subject of 
renormalization in the presence of evanescent operators in Section 4. In this section we present 
results for the bare matching coefficients. 

We define the Fourier transform of the coefficient of the SCET n operators O r ,i(s,t) as 

DtM U)= jds dt e -^vt e iusn.P t) ^ (lg) 

where L>r t i(s,t) is the position-space coefficient. The jet functions relating bare coefficients 
are defined via 

= ^£dvJi(u,v,g 2 (2E U )- e )C^ n (v), (19) 

which we expand as 

J t (u,v,g 2 (2ELu)-*) = J^f -SaSiu-v^S + O) 
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Here g is the bare coupling; the dependence of J± on g, E and a; is determined by dimensional 
analysis and invariance under rescaling of the light-cone vector n. The symbols S and O 
represent the color-singlet and color-octet configurations of four-quark operators. Restoring 
the color indices which are to be contracted with the operators in (15), we define: 
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where for convenience we include the tree-level factors in the definitions of S and O. 

The one-loop jet-functions for the three operators in (17) are obtained from the SCETi 
diagrams shown in Figure 3. Diagrams involving soft modes need not be considered, since 
these momentum regions are reproduced exactly by the corresponding SCETn contributions, 
and hence cancel in the matching. We proceed by evaluating the loop diagrams for an arbi- 
trary coefficient function Cr in (14). To perform the integration over the ci-dimensional loop 
momentum k, we first integrate over the n ■ k component using the method of residues, and 
then perform the k± integral for the d — 2 perpendicular dimensions. The remaining integra- 
tion over n ■ k takes the form of equation (19), from which we read off the contributions to the 
jet-functions. In this way, we find 
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(24) 



The plus distribution appearing in is defined for symmetric functions g(u,v) to act on 
test functions f(v) as 



J dv [g(u, v)]+ f(v) = J dvg(u, v) [f(v) - /(«)] . 



(25) 



We now return to the particular cases of the scalar, vector and tensor current operators 
and extract their jet-functions from the general expressions above. Through one-loop order, 
one physical and one evanescent four-quark operator are sufficient for the scalar case: 



S i = 7± ® l±a , 



S 2 = if 1 ® l±af3y ■ 



(26) 
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To the same order, there are four physical and four evanescent vector operators 
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and a grand total of fourteen tensor operators 
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(28) 



of which the last eight are evanescent; in particular, Ot7, though appearing already at tree- 
level, is evanescent. 2 As is evident from the above equations many of these operators are 
directly related to each other, e.g. Ovi to Osi, etc., and the matching onto these operators 
involves the same jet-functions. To obtain the jet-function in a particular case, we first rewrite 
the SCETj current in the general form (13), and then evaluate the three corresponding op- 
erators of (17). Matching to the basis of four-quark operators defined in (26), (27) and (28) 
yields the jet-functions in terms of the functions J7i defined in (19). The independent matching 
relations are obtained from those operators with no factors of v M or n M , and with zero, one or 
two perpendicular Lorentz indices: 
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(29) 



2 Although three independent structures are possible for strings of five 7j_ matrices with the appropriate 
antisymmctrizations, only Otis and Oti4 arise in the matching. 
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Since they are related to the operators above simply by overall factors of and n**, the 
remaining matching relations onto CVi,2,5,6 an d onto Oti,...,5,8,...,i2 may be obtained from (29), 
using (26), (27) and (28), and the analogous relations for the SCETj currents in (7) and (8). 

From (29), we note that for a particular SCETi current Jr in (13), the index contrac- 
tions in the resulting SCETn operators Or, 2 and 0r,3 in (17) lead to e-dependent prefactors, 
which need to be expanded around d = 4 before renormalization in order to arrive at MS 
subtracted jet-functions. Before performing this renormalization, we first consider carefully 
the renormalization scheme dependence implicit in our above choice of evanescent operator 
basis. 

4 Evanescent Operators in SCET 

We have seen in the previous section that in dimensional regularization, an infinite tower of 
evanescent operators is generated upon matching SCETj onto SCET n . These operators can 
mix into the physical basis, and so their effects must be included for a consistent analysis. To 
ensure that these operators do not introduce additional nonperturbative parameters into the 
theory, it is important to show that a renormalization scheme can be chosen such that their 
renormalized matrix elements vanish. Here we show that this is in fact the case. A related 
issue concerns the choice of evanescent operator basis. As we will discuss, the renormalized 
matrix elements of physical operators depend on this basis choice, so that choosing a particular 
basis of evanescent operators induces a renormalization scheme dependence in the physical 
operators. For the four-quark SCETn operators defining the hard-scattering part of the form 
factors, we choose the basis such that the resulting renormalization scheme corresponds to 
the MS scheme after Fierz transformation. In Section 4.1 we present a general discussion of 
evanescent operators in SCET, while in Section 4.2 we apply the general formalism to the 
operators appearing in the form factor analysis. 

4.1 Evanescent operators and the Sudakov problem 

Evanescent operators first appeared in the case of local four-fermion interactions [18]. At 
one-loop order, it is easy to see that the Green functions of bare evanescent operators are 
finite and local, and hence can be made to vanish by the addition of a suitable counterterm. 
The finiteness of the Green functions results from the 1/e pole from the momentum integral of 
the corresponding Feynman diagrams being cancelled by a factor e in the Dirac algebra; this 
latter factor of e always appears when evanescent operators "cross the threshold" and mix into 
physical operators, reflecting the fact that evanescents, and hence any contractions involving 
them, vanish in d — 4. The locality of these particular Green functions is a consequence of 
the locality of the coefficient of the 1/e pole of one-loop integrals. A renormalization scheme 
in which matrix elements of evanescent operators vanish to all orders in perturbation theory 
may then be established by induction [19]. In the present situation, where Green functions of 
the (nonlocal) SCET current operators contain 1/e 2 double poles already at one-loop order, 
one may ask if an analogous renormalization scheme for evanescent operators exists. We will 
show that this is indeed the case, both for SCETj and SCET n . In practical applications where 
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we will ultimately match onto SCET n , we may simply use the conventional MS scheme also 
for evanescent SCETt operators and perform finite subtractions only in the final low energy 
theory. This is the procedure adopted in Section 4.2. 

Let us consider first the simpler case of the evanescent SCETi current operator, J^ 5 , 
denoted by Jt5 in the following. Under renormalization, this operator mixes with J^, denoted 
by Jti- 3 We proceed to isolate the linear combination of J T \ and J T5 which corresponds to 
the physical renormalized operator. It will be convenient to define the following notation for 
these operators, considered as bare quantities, 

di(s,r) = J^ v (s,r) = X(sn) A ±a (rn) ^ u h(0) , 

a 2 ( s , r) = X(sn) 4 ± {rn) 7^7^ h(0) = -2e X(sn) A ± (rn) 7 f h(0) = -2e J^^(s, r) . (30) 

From the Feynman rules of SCETt and the projection properties of the spinor field X, it follows 
that the operators Si and 82 close under renormalization [16] . Defining renormalized operators 
via dT n = J2j Zij Jj are , the anomalous dimension matrix takes the form 

(31) 

where, in the MS scheme, it is understood that only pole terms in the dimensional regulator 
e = 2 — d/2 are kept in Z — l. 4 If we now consider a Green function with an insertion of 
the renormalized operator, ^ en , the result will be the sum of two terms, the first proportional 
to the tree-level matrix element of Si, and the second proportional to the tree-level matrix 
element of Si- The renormalization constants of Z ensure that the coefficients for both terms 
are finite. But Si is evanescent, and S2 contains an explicit factor of e, so their tree-level matrix 
elements vanish at d = 4; since the coefficients have been made finite by renormalization, the 
final result also vanishes in d — 4. Returning to our original operators, we define 




Jt7 ) = ( 3[ cn \ = Z u -2eZ l2 
J tiJ ~ \-TedTj V 222 ) \ J tT ) 




(32) 



J™i and Jy 5 n are equal to their bare counterparts at tree level, and have finite Green functions 
to all orders in the loop expansion, so that (32) defines a valid renormalization scheme. Clearly 
J|? 5 n , being equal to S r \ n , has vanishing matrix elements in this scheme. We may note that the 
renormalization of the new basis involves a finite subtraction from the term eZ\2, even if the 
original matrix Z in (31) is defined in the MS scheme. 

None of this argument depended on the order in the loop expansion, or on the choice of 
renormalization scheme used to renormalize Si and Si- Also, since the renormalized matrix 
element of Jt5 vanishes regardless of renormalization scale, we find that 

= — (Jt5) = -1t5,T5(Jt5) — 1t5,ti(Jti) = — 7t5,ti(-/ti) , (33) 

aln/i 



3 This mixing only begins at two-loop order. 

4 Z 2 2 is given by the pole part of Z\\ + 2(1 — e)Z\2 [16]; this fact will not be relevant to our analysis. 
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Figure 4: Abelian contribution to the anomalous dimension of the 5-type SCETi current 
operators. Full lines denote soft fields, dashed lines hard-collinear fields. 

where the last equality follows from the vanishing of (Jts). Since the matrix element of 
the physical operator is nonzero, it follows that 7t5,ti = 0. The evolution equation for the 
coefficient function of the physical operator is therefore simply 

— Cti = 7ti,tiCti • (34) 

In particular, the evolution equation is independent of CV 5 , and the coefficient of the renor- 
malized evanescent operator is irrelevant in the calculation of physical processes. 

The simplicity of these arguments hides some interesting implications. For instance, the 
vanishing of 7t5,ti implies that at one-loop order eZyi is finite, or that Zyi has at most a single 
pole at e = 0. The one-loop 1/e 2 double poles must occur only in the diagonal elements of Z, 
which do not mix the physical and evanescent operators. To illustrate how this works, we recall 
that the anomalous dimensions of the -B-type SCETi current operators may be extracted from 
the UV divergences of 3-point Green functions. For simplicity we consider the abelian case, for 
which the relevant diagrams are shown in Figure 4. The soft contribution is represented by the 
first diagram, while the remaining diagrams give the hard-collinear contribution. Regulating 
infrared divergences by evaluating the Green function at an off-shell value of the hard-collinear 
quark momentum p, we find 

softocf^y" (-\+.. 

hard-collinear oc f — ) ( -r + ■ ■ • } . (35) 

By general renormalization arguments [24], the coefficient of the 1/e pole in the one-loop 
3-point QCD Green function, which the effective theory is representing, is local and hence 
independent of p 2 . In fact, in SCETi, the QCD result is reproduced by including the hard 
contribution, accounted for by the Wilson coefficient, in addition to the hard-collinear and 
soft contributions. This hard contribution is independent of p 2 , but may depend on the scales 
m b and E, so that the SCETi divergences may be nonlocal over distances ~ 1/E. This is 
exactly the structure obtained after taking the sum of the two contributions in (35): 

1 2 *2E 

soft + hard-collinear oc — In h • • • , (36) 

e 2 e /i 

where the remaining divergences are local at the level of SCETi, i.e. independent of p 2 . Since 
the Dirac structure of the original current is unaffected by the soft gluon exchange, the 1/e 2 
term in (36) contributes only to the diagonal terms in Z. 
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We may present similar arguments to show that the renormalized matrix elements of the 
evanescent SCETn operators of the form (15) vanish in d — 4. Since the heavy quark couples 
at leading power only through a spin-independent v ■ A s interaction, no new Dirac structures 
can be induced by the soft sector. Similarly, the soft-collinear modes couple to collinear 
modes only through n ■ A sc , and to soft modes only through n ■ A sc . All evanescent operators 
therefore arise from the collinear sector, and we will only write the collinear part of the four- 
quark operators in our discussion. For simplicity, we restrict attention to the color-singlet 
case, 5 and study the operators, 



[OiU= [Xf ± ] a [T ± X] 



b ' 



(37) 



with spinor indices a, b. Here T±, T±, represent as-yet unspecified Dirac structures (T± is not 
related to Tj_). It follows from the Feynman rules of SCETn that 0\ mixes with the infinite 
tower of operators, 



[o n u= [x 7 r 



'7/ 



71 — 1) 



• DC 



j b 



(38) 



with O n appearing at n-th order in perturbation theory. Here the indices z/j are not anti- 
symmetrized. By general renormalization arguments [24], renormalization constants Z nm may 
be chosen so that with 0™ n = ^ m Z nm 0^ re , Green functions of the renormalized operators, 
0^ en , are finite. 

In particular, in the representation of the scalar QCD current, we may consider the first 
evanescent operator, Os2 in (26), for which we take f j_ <g> T± = •y^ pa <g> 7± Mp(7 . If we also 
introduce the operator appearing at tree-level, 



[o ] ab =[x^ ± ] a [ 7±M x] 6 , 

then the renormalization matrix takes the form, 6 



(39) 



/ O T cn \ 

rcn 
Qren 

V : I 



(Z 



00 ^01 ^02 
Zu Z12 

Z21 Z22 



\ ( oY c \ 



\ ■ 



bare 



/Obarc 
w 2 



(40) 



The reasoning now proceeds as in the previous case. The operators O n for n > 1 close under 
renormalization, and after renormalization, a Green function involving the insertion of one of 
these operators is a linear combination of finite coefficients times tree-level matrix elements. 
Since the tree-level matrix elements vanish in d — 4, the renormalized Green functions of the 
evanescent operators 0^ en , n > 1, also vanish. Again, since this is true for any renormalization 

5 For the color-singlet case, the soft-collinear modes decouple entirely [23], so that the collinear and soft 
sectors factorize. 

In fact, since at higher loop order the number of factors of 7^ may not decrease, the matrix of renormal- 



ization constants is block triangular, with Z n 
discussion. 



for n < m. We do not make use of this fact in the present 
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point, the elements of the anomalous dimension matrix describing the mixing of the physical 
operator into evanescents vanish. Although we have specialized here to the case of the scalar 
current, the arguments may be easily extended to the general case. 

The basis of evanescent operators O n is of course not unique. We could for instance have 
chosen the antisymmetrized basis, 



[0„U= [X 7 f •••7l 2n+ll ] a [7±^---7^ + A, 
or any other basis which is related to the original basis of operators by 



(41) 



Oo = Q 



and 



O n =^ A nrnO m + eB n 6 for U > 1 , 



m=l 



(42) 



with coefficients A nm , B n , finite at d — 4. The same renormalized operators may be expressed 
in terms of the new basis of bare operators as 



/ O T cn \ 

rcn 

oi en 



( Zqo + e Em=l ^Om-Dm 



m=l 

E°° 7 4 
m= l ^lm^ml 

C Em=l Z 2mB m Y2m=l ^2m^ml 



Sm=2 ^0m^m2 



■VOO 



E m =2 Zi m A m2 ■ ■ • 
E^=2 ^2m^m2 * " " 



V ; / V 



^bare 
Qbare 



(43) 



■J \ J 



If we identify the resulting matrix with the renormalization constant matrix Z for the new 
basis, then we find for n > 1 the typical counterterms, Z n o = ^Ylm=i Z n mB m , involving finite 
constants which must be included in order to ensure the vanishing of the renormalized evanes- 
cent operators. We also find, however, a finite contribution to Z o = Z 00 + e Ylm=i z omB m - If, 
as is customary, we impose the MS scheme for the physical operators and allow finite subtrac- 
tions only for evanescent operators, in the terms Z n0 for n > 1, then the finite contribution to 
Z Q o must be absorbed into the Wilson coefficient function of Oq CI1 . Since this finite contribu- 
tion depends on the choice of basis in (42), a scheme dependence is induced in the coefficient 
functions of the operators; this scheme dependence cancels against the corresponding scheme 
dependence in the operator matrix elements. 



4.2 Scheme dependence of the jet functions 

From the preceding discussion, we find that using the MS scheme for the physical operators, 
a change of evanescent operator basis corresponds to a change of renormalization scheme. 
In this section we work in the opposite direction, and ask if it is possible to find a basis of 
evanescent operators which yields a given renormalization scheme for the physical operators. 
In particular, when we take the matrix elements of physical operators describing form factors, 
we would like the LCDAs appearing in (1) to be defined in the MS scheme. 

We start by defining a Fierz rearranged operator 0[ for each physical operator Oj in such 
a way that the bare operators coincide in four dimensions: 

S i = XilK Q 1±a X <- 0' S1 = 2 X L |x Q L |lK + 2 X R | X qJ^K , 
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Figure 5: One-loop graphs contributing to the difference between the matrix element of the 
original and the Fierz rearranged four-quark operators. 

V 3 = X 7 ^7l^ Q7± a X <- 0' V3 = 2X L7 ^|x qJ^H + (L <- R) , 
Ova = tilipi Q 1±a X <-> 0' V4 = 2X L |x Q^l^K + (L ^ R) , 
T6 = XHIZ'K Q 7 ± a X <- 0' T& = 2X L |x Ql^'K + (L <- R) , 

T 7 = X 7 7^ Q 7±Q X <- 0^ = -2X L7 i M |xQ fi | 7 ^+(L^i?). (44) 

The remaining physical vector and tensor operators are obtained by multiplying the above 
operators by factors of and n M . In fact, since Ova and Ot% are obtained from Osi by 
insertion of Dirac structures next to the heavy quark field, these operators don't require an 
independent analysis. Similarly, from the identity — 2 7 " Miy = 7 x7±7± — 7l7±7j\ we niay relate 
Ot7 to Ov3- It is convenient to treat Ot7 in the same manner as the other operators in 
(44), even though it is evanescent. We will find that the resulting renormalization scheme 
automatically ensures the vanishing of the renormalized Ot7- 

We define the renormalized physical operators Oi and 0[ in the MS scheme. The heavy-to- 
light meson matrix elements of the Fierz rearranged operators 0\ will then yield the product of 
the LCDAs of the two mesons in the MS scheme. Our task is to define the evanescent operators 
in such a way that physical matrix elements are the same for Oi and 0\. For instance, if we 
consider a four-point Green function involving 0' SI , then at one-loop there are contributions 
proportional to the tree-level structure, as well as a contribution proportional to 

\ x = \ x (4 - 86 + 4e 2 )| - - 8 | . (45) 

This contribution arises from the second diagram in Figure 5. In the MS scheme, only the 
1/e pole is subtracted to obtain the renormalized result. Similarly, the matrix element of 0$\ 
depicted in the first diagram in Figure 5 contains the same contribution proportional to the 
tree- level structure 7" <g> 7 ^ Q , and a contribution proportional to 

~ e x 7^7±7l ® 7± 7 7±/37±a = - x (4 - 6e)7^ <g> 7 ±a - ^ x if 1 <g> 7±Q/37 

^-6 7 !®7± a - (46) 
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Using MS subtraction in the original basis (26) yields a renormalized result different from the 
corresponding result in (45). In place of (46), we should instead write 

~ e x 7^7?7l ® 7± 7 7±/37±a 

= \ x (4 - 8e + 4e 2 ) 7 ^ ® 7±Q - 1 x [ 7 fr <g> 7±Q/?7 - 2e(l - 2e) 7 ? <g> 7 ±J 

^-8 7 !®7± a - (47) 

The resulting contributions to the renormalized matrix elements of 0' SI and Osi will therefore 
agree if we choose in place of Os2 in (26), 

S 2 = S2 - 2e(l - 2e)0 51 . (48) 

Beyond one-loop, the other evanescent operators Oss, Osa, . . . are likewise determined. Pro- 
ceeding in the same way for the vector case, we find 

Ovi = V7 + 2e(l + 2e)0 V3 - 4eO V4 , 

V8 = O y8 -2e(l-2e)Oy 4 . (49) 

Dropping the C(e 2 ) terms which are irrelevant through one-loop order, the scalar and vector 
jet-functions in this new basis are (cf. (29)) 

Js 

Pj. Osi (Ji + 2(l-e)j 2 + 2ej 3 \ (50) 

JS Os2 \ JS 



jB' 
J VI 



jB> 
J VA 



Jv = 



o V3 

d V7 
o V8 



( Ji~2eJ 3 \ 

J 1 + 2(l-e)J 2 + 2eJ 3 2eJ 2 + 4eJ 3 

Js 

Js 



(51) 



Although they will not be relevant to the form factor analysis, for completeness we give here 
also the new basis of evanescent operators for the tensor case, 



Otis = T13 - 2e(l - 2e)0 T6 , 
Otia = T1A + 2e(l + 2e)0 T7 + 4eO T6 , 
and the corresponding jet functions, 



(52) 



jB> 



jB' 



J T 



Ot6 

Otis 
Otia 



( J x + 2(1 - e)J 2 + 2eJ 3 -2eJ 2 - 4eJ 3 \ 





Js 




Ji ~ 2eJ 3 


Js 



(53) 
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We may recall from (44) that the operators Oti and 0' T7 vanish at tree level in d = 4. In 
fact, since it is multiplicatively renormalized, the renormalized 0' T7 then vanishes to all orders. 
Using the new basis of evanescents, Oti3, Ari4, • • • , the same is true for Ot7, since the new 
basis was chosen precisely to ensure that renormalized matrix elements of O77 and 0' T7 agree. 
In this scheme, therefore, Ore is identified with the renormalized physical operator, and the 
renormalized Ot7 vanishes. It is interesting to note that the evanescent J|? 5 gives a finite 
contribution at one-loop order to the bare matching coefficient of the physical operator Ox% 
in (53). As discussed in Section 4.1, since we use the MS scheme for all operators in SCETi, 
the corresponding matching coefficient between renormalized operators will also be nonzero. 

Let us illustrate the above arguments by an explicit evaluation of the diagrams in Fig- 
ure 5 for the scalar case. In terms of tree-level matrix elements, the difference between the 
renormalized matrix elements at one-loop order is 

(O r s7(u))(u') - (0' s T(u))(u') = J dvA(v,u>) [(4 - 6e)(O sl (u))(v) - (0 S2 (u))(v) 

-A(l-e) 2 (0' sl (u))(v)}+ tdvZ sl>S2 {u,v){0 S2 {v)){u'). (54) 



Here u is the Fourier transform variable in (18) and u' is the momentum fraction of the outgoing 
collinear quark in the diagrams in Figure 5. At tree level, we have (Oi(u))(u') oc 5(u — u'). 
The function A is given by 



A(u,u') 



a s Cp 1 
4tt 2e 



u 



u 



-6{u' -u) + —6(u - v! 



u 



W 



(~U) T (l)f, 



An u 



^j^K, M ) + o( e °), 



(55) 



where in the last equality we compare to the explicit form of in (24). Since the divergences 
have to cancel, we see that A = Z$i t s2, and in our original evanescent operator basis the finite 
difference in the renormalized matrix elements becomes 

(O r s7(u))(u') - (0'™(u))(u') = J dv2eA(v,u')(O sl (u))(v) . (56) 
The one-loop Wilson coefficient functions of these operators must therefore fulfill the relation 
J^(u,v) - Jsi(u,v) = j dyJlT(y,v)2eA(y,u) = 2eJ 3 (u,v) , (57) 

which is in agreement with (29) and (50), since Jsi = J' sv 

Having determined the operator basis, we now renormalize the jet-functions. After MS 
subtractions, the quantities J\ + 2(1 — e)J 2 + 2eJ" 3 and J\ — 2ej7" 3 yield J\\ and J±_ from [16], 
respectively. We recall that J\\ is the jet function appearing in all form factors describing B- 
decays into pseudoscalar mesons, or longitudinally-polarized vector mesons, while J7j_ appears 
for all perpendicularly-polarized vector mesons. In [16] we introduced the following notation, 



J\\,±.{u,v,L,n) 



47rCi?Q; s (/i) 1 



N 



2Eu 



-8{u 



47T 



j l± (u,v,L) + 0(a 2 s ) 



(58) 
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where L = \n(2Eu / ji 2 ) , jy = j 1 + and jj_ = ji — J3. In terms of the one-loop expressions 
(22) and (23), after MS subtractions we have 



(i) 



2(1 -eU 



(i) 



2eJ. 



(i) 



Ji 

J2 
J3 



(59) 



The quantity 2ej7~2 also appears, and we find 



2e4 1 \u,v) 



—6(u -v) + —9{y - u) 
uv uu 



uv — (1 - u — v) 2 . _ , 

H ^— —0(1 -u-v) 

uuv 



+ 2C, 



v-u vu 

— —0(v — u) 

uu u 



• (60) 



However, this jet function occurs only in the matching onto operators irrelevant to the form 
factor analysis. Similarly, the quantity J% appears in the matching onto evanescents. 



5 Application to heavy-to-light form factors 

In this section we study the numerical size of the matching corrections in the factorization 
formula for the form factors. To begin, we review the situation for the A-type currents. 
This case is especially simple because no convolution over momentum fractions is involved, so 
that no assumptions on the shape of hadronic wave functions are necessary in order to make 
statements about the size of one-loop corrections relative to the tree-level values. We consider 
the combinations of Wilson coefficients for the A-type currents which yield the Cj in the first 
term of (1) for the form factors F + , F , F T) V and T\. To distinguish these from the Cf in the 
hard-scattering terms, we will denote them by Cf in the following. Form factor normalizations 
are chosen such that Cf is equal to unity at tree level. The coefficients for the form factors 
of longitudinally-polarized vector mesons (A , V — A 2 , T 2 — (2E/mg)r 3 ) are identical to the 
pseudoscalar case. Also, for perpendicularly-polarized vector mesons, form factors V and T 2 
are related to A 1 and Ti respectively. The same relations also hold for the hard-scattering 
terms. Using the expressions from [5, 17], we find that the one-loop corrections are remarkably 
similar for all form factors. As illustrated in the plot in Figure 5, the corrections are negative, 
ranging in size from approximately 10 — 15% at E — mi>/2 to 20 — 25% at E — m^/A. The 
symmetry breaking part of these corrections is very small, less than 5% in the energy range 
considered. This may indicate that the bulk of the one-loop correction is simply due to the 
choice of renormalize scheme or factorization scale. 

The analysis is more complicated for the 5-type current operators. Since the hard- 
scattering kernel T^E, ou,u) appearing in the factorization formula (1) depends on variables 
lo and u, the numerical impact of the loop corrections can be assessed only after convolution 
over these variables. Knowledge of the wavefunctions of the 5-meson and light final state 
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Figure 6: One-loop coefficients Cf(E) of the soft-overlap terms in the heavy-to-light form 
factors as a function of the recoil energy. The QCD as well as the effective theory renormal- 
ization scale are set equal to mj, and we take a s (mb) = 0.22. All coefficients are equal to one 
at tree level. 

meson appearing in the convolution is therefore necessary. To begin, we recall the form of the 
hard-scattering term [16] 



where Jr = J\\ for pseudoscalar and longitudinally-polarized vector mesons and Jr = J± for 
perpendicularly polarized vector mesons. The quantity K F (n) = 1 + ^^(3m^ - 2) relates 
the QCD and the HQET fi-meson decay constants. In order to avoid large perturbative 
logarithms in the coefficients Cj 3 , one would like to set /i 2 ~ m 2 , whereas the scale associated 
with the jet-function is lower, \i ~ ^/rribAh, where Ah is a typical hadronic scale. Since the form 
factor involves different scales, we see that we cannot avoid large perturbative logarithms in 
a fixed order calculation. These logarithms can be resummed by solving the renormalization 
group equations for the Wilson coefficients. This is discussed in detail in [16], where we 
resummed all leading logarithms of the different energy scales in the problem. The two effects 
we wish to study here are the one-loop corrections to Cf , and to Jr- We study the two 
corrections in turn, evaluating each of them at their natural scale /x, but remind the reader 
that such a scale choice is strictly only possible after renormalization group improvement. 
Let us start by considering the hard-scale matching corrections to the B-type Wilson 
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coefficients. The five independent combinations appearing in the form factors are 



— (_/y 2 -+- (_/y 3 — 1 h . . . , 

Cf t = = 1 + • • • > 

C B _ r B> _ n , 

u Ai — °y4 — u + • • ■ 5 



At tree level they take the values indicated; the pertinent one-loop expressions were given in 
Section 2. To evaluate the impact of these corrections on the form factors we work with the 
tree-level jet-function, which is independent of uj: 

rr i \ 4irC F a s ((j,) 1 2 

Jt{u,v) = — <J(u - u) + ©(aj . (63) 

The size of the hard-matching corrections is then determined by convoluting the Wilson co- 
efficient Cf with the light-meson wave function. Using for simplicity the asymptotic form 
</>m(w) = 6mm of the light meson LCDA, the relevant integral is 

C? (cff) (E, W ) = J\uu Cf(E, u, W ) . (64) 

These effective coefficients coincide with the Cf at tree level, since dependence on the mo- 
mentum fraction u appears only at one-loop order. In terms of the correction to the 
form-factor is 



— -m 



3irC F a s (n) f B f M (v) 



Cf (cff) (£,/,), (65) 



N m B X B (jj,)_ 

where As(/i) is defined as the first inverse moment of the 5-meson LCDA, 

A bV) = / — <M^,/i). (66) 
Jo w 

The quantity inside square brackets in (65) coincides, up to RG factors, with the tree-level 
value of the function Hm introduced in [16]. In Figure 7, we plot the one-loop contributions 
to Cf^ cS \E, «). Choosing the scale « = m b , we find that the one-loop corrections are of order 
20% of the tree- level values given in (62). 

Finally, to isolate the one-loop corrections from the jet functions, we use tree-level matching 
for the hard coefficients. Since these tree-level coefficients are independent of the momentum 
fraction u, in this approximation the resulting convolution integrals in (61) are universal to 



24 



0.3 




0.2 



0.3 



0.6 



0.7 



0.8 



0.9 



1 




Figure 7: One- loop corrections to the form factors from the Wilson coefficients Cf(E, u). The 
figure shows the one-loop correction to the quantity (E) that arises after convolution of 

the coefficient Cf(E,u) with the leading order jet-function and the asymptotic pion LCDA, 
see (64). The renormalization scale is set equal to m b in both QCD and the effective theory. 

all form factors involving the same final-state meson. For pseudoscalar and longitudinally- 
polarized vector mesons, the result is proportional to [16] 



where for simplicity we have again used the asymptotic form for the light-meson LCDAs. We 
take Uf = 4 as the number of light quark flavors. For perpendicularly-polarized vector mesons, 
where J\\ is replaced by J7j_, the coefficient of the double logarithm is the same as in (67), the 
coefficient of the single logarithm is changed from (—19/3 + tt 2 /9) to (—6 + 7r 2 /9), and the 
constant term is changed from 3.99 to 0.73. Here the angle brackets denote averages over the 
5-meson LCDA, 



Jo w 

To investigate the size of the one-loop corrections, we may take illustrative models for the 
£>-meson LCDA. For example, in [25], 




(67) 




(68) 



4>b{w,Hg) 



A B 



\ B = - (m B - m b ) w 0.32 GeV , 



(69) 



25 



where /x is a low hadronic scale at which the model is assumed valid. For the averages in 
(67), we evaluate the wavefunction appearing in (68) at the scale /j,q, and take the intermediate 
scale \x\ = rribA h , with A h = 0.5 GeV, as the argument of the jet function. For maximal energy 
E = mi>/ 2 this model yields 



K / V K 
/ 2Euj\\ 2 /'2-EAsexp 



7E\ n 2 



ln \-w))^{-^^) + ^- 2 - 7 - (70) 

The resulting hard-scattering corrections in (67) are approximately +30% and +40% for J\\ and 
J7j_ respectively. Neglecting the implicit model-dependence of these results, the uncertainty 
due to nib and hence of is still significant; for instance, taking \b = 0.46 GeV [26] reduces 
the size of corrections to approximately +20% and +30%. Using the model wavefunction from 
a more recent analysis in [26], with the same value of the intermediate scale in the jet function, 
and at maximal energy E = nib/2, gives 

In ^ ~ -0.7 ± 0.4 , ^ln 2 ^ ~ 1.0 ± 0.6 , (71) 

where the wavefunction is evaluated at the scale 1 GeV. In this model, the corrections for J\\ 
and J7j_ are approximately 15% and 25% at E — nib/2, and grow larger at small values of E, 
reaching 30% and 40% at E = m b /A. In evaluating the averages in (67) with wavefunctions 
evaluated at the low scale /i ~ 1 GeV, and at the same time taking /ij = m b A h as the argument 
of the jet functions, we have neglected the RG factors relating the scales /x and /ij. From [16], 
we find the resulting correction to (67) for this choice of scales to be approximately —15%. 



6 Discussion and conclusions 

We have presented in this paper the complete results for the one-loop matching corrections to 
the hard-scattering coefficients in the factorization formula (1) which describes heavy-to-light 
form factors at large recoil. Our results (9) and (10) for the hard-scale matching coefficients 
Cf in the vector and tensor case extend our previous scalar current results reported in [16], 
and are in agreement with an independent calculation [17]. Results for the intermediate-scale 
coefficients, J7j| and J7j_ in (58), which are relevant to B decay form factors were reported in 
[16]; in this paper we have presented general results, and details of the calculation. 

An interesting theoretical issue involves the evanescent operators which arise upon match- 
ing QCD onto the effective theory. We showed that even in the presence of divergences of the 
Sudakov type, a renormalization scheme exists for which renormalized evanescent operators 
vanish. The issue arises already in SCETi, for the -B-type operators, and we discussed this 
case first. This example is especially interesting since the mixing involves only two operators, 
in contrast to the usual infinite tower of evanescent operators induced by radiative correc- 
tions to four-fermion interactions. The case of the nonlocal four-quark operators of SCETn 
was then treated along the same lines. On a practical level, use of such a renormalization 
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scheme implies that matching coefficients onto evanescent operators become irrelevant below 
the matching scale, since the renormalized operator matrix elements vanish, and the physical 
operators do not mix into the evanescents. We also discussed the choice of evanescent oper- 
ator basis. Different bases correspond to different renormalization schemes for the physical 
operators. We isolated a particular basis which corresponds to the MS scheme after Fierz 
transformation. 

We have presented an analysis of the phenomenological impact of one-loop radiative cor- 
rections in Section 5. Using the known results for the A-type currents [5, 17], we find that the 
matching corrections to the soft-overlap terms are remarkably universal to all form factors. 
The corrections are of order 15%, but differ in all cases by less than 5%. This effect seems not 
to have been noticed previously in the literature, and implies that any breaking of the naive 
large-energy spin-symmetry relations by more than 5% is due to either the hard-scattering 
contributions, or to uncalculated power corrections. We have also investigated the numeri- 
cal impact of radiative corrections to the hard-scattering coefficients, finding the hard-scale 
corrections to be < 20% for all form factors. These corrections determine the size of viola- 
tions to the universality of the hard-scattering terms for different form factors which holds 
at tree level for the hard-scale coefficient. Radiative corrections to the jet functions are less 
important from the point of view of universality, since they contribute identically to all form 
factors involving the same final state meson. On the other hand, for the purpose of calculat- 
ing the hard-scattering contributions in terms of the B meson wavefunction, the jet function 
corrections are expected to be more important than the hard-scale corrections, since the cou- 
pling constant a s (fi) is larger at the intermediate scale jj 2 ~ Aqcd^- For this reason it has 
been speculated that the jet-functions might not have a convergent perturbative expansion 
for physical values of m b [27]. We evaluated the perturbative corrections to the jet-functions 
for two models of 0s (tc), finding that one- loop corrections are of order 20 — 30%, with a large 
uncertainty owing to the uncertainty in the model parameters. 

We have focused here on the heavy-to-light meson form factors, but the same Wilson 
coefficients are also relevant for other processes. The jet-functions calculated here appear in 
other exclusive processes, such as B — > K* r y and the decays to two light mesons [27]. The 
one-loop matching corrections to the SCETi current operators might become relevant for the 
study of power corrections to inclusive decay distributions for processes like B — > X u tv and 
B^X sl [28]. 
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